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Рассмотрены проблемы нелинейной динамики процесса об- 
работки материалов резанием. Предложена математическая 
модель динамической системы, учитывающая динамическую 
связь, формируемую процессом резания. Динамическая под- 
система инструмента представлена линейной динамической 
системой в плоскости, нормальной к поверхности резания. 
На основе выполненных исследований даются рекомендации 
по проектированию систем, имеющих требуемые стационар- 
ные многообразия в пространстве состояния. Полученные 
данные позволяют сделать следующие выводы. Динамиче- 
ская система резания может служить примером сложного 
поведения, зависящего от параметров подсистемы инстру- 
мента и представления сил в координатах состояния. Пред- 
ставление сил в координатах состояния зависит от множества 
физических процессов, сопровождающих обработку. Важно 
отметить, что в зависимости от моделей, взаимодействующих 
через процесс обработки, в окрестности равновесия могут 
формироваться различные критерии множества, вызывающие 
регулярные или нерегулярные свойства формируемой при ре- 
зании поверхности. 


Ключевые слова: процесс резания материалов, динамиче- 
ская система, инвариантные многообразия, бифуркации. 
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Введение. Проблемы динамики процесса резания исследуются в течение последних 50 лет. При этом внимание уделя- 
лось, главным образом, двум вопросам: условиям и механизмам возбуждения автоколебаний [1—5] и анализу устойчи- 
вости процесса резания [6-10]. Довольно активно изучаются подсистемы инструмента и обрабатываемой заготовки, 
их взаимодействия через динамическую связь, формируемую процессом резания. Полученные в результате данные 
служат основой для исследования динамики процесса резания. 

Для раскрытия основных свойств динамической системы можно использовать базовую динамическую модель 
процесса обработки, которая рассматривает упругие деформационные смещения инструмента в плоскости, нормаль- 
ной к поверхности резания. При этом динамическая связь, формируемая процессом резания, представляет модель сил 
резания в координатах состояния системы [11-15]. 
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В данном случае учитываются: 

— нелинейная зависимость сил резания от технологических режимов [16-18]; 

— периодические изменения параметров динамической системы [ 19—21]; 

— закономерности формирования новых поверхностей [22—25]. 

При рассмотрении устойчивости системы с периодически изменяющимися параметрами, а также при форми- 
ровании новых поверхностей резания используется теория Флоке [22]. 

В статье развиваются указанные выше представления о динамической системе резания. Основное внимание 
уделяется вопросам динамической самоорганизации системы, а также ее бифуркациям в параметрическом простран- 
стве. При этом не принимаются во внимание изгибные деформационные смещения инструмента. 

Базовая математическая модель. Если не принимать во внимание зависимости сил от координат состояния, то 
уравнение динамики процесса можно представить в виде [14]. 
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где т = И В ое а. ” | — соответственно матрицы инерционных, диссипативных и упру- 
О т Й, › й, › С2  С22 
гих коэффициентов подсистемы инструмента; Х ={Х.,Х | — вектор упругих деформационных смещений инстру- 
р) — 1? 2 у руг р тру 


мента в плоскости, нормальной к поверхности резания; Ё(!) ={А (1), Е, (#) }’ — вектор сил резания. 
Упругие деформационные смещения рассматриваются в системе координат, двигающейся вдоль заготовки со 


скоростью движения суппорта \У.. Угловая скорость вращения заготовки считается постоянной. Она определяет по- 
стоянную составляющую скорости резания У. Поэтому в установившемся состоянии величина подачи на оборот 
ь ‚ = СОП51. 

Сила резания формируется в виде суммы двух составляющих, действующих на переднюю грань инструмента 
Е‘ и на его заднюю грань Ё (2) ‚то есть Е=Е® + Е (рис. 1). При заданной геометрии инструмента технологиче- 


ские режимы (глубина резания {, и величина подачи на оборот 5, ) однозначно определяют ширину срезаемого слоя 
Ь и его толщину @,т. к. а=5ьзшф,ар= Г, /зшф ( Ф(рис. 1) — главный угол инструмента в плане). Деформа- 


ционные смещения в направлении, нормальном к плоскости (Х, — Х,) ‚ Не принимаются во внимание, Так как де- 


формации В этом направлении практически не вызывают изменения площади срезаемого слоя № : Кроме того, при 


точении обычно выполняется соотношение 5 ь ((1 р: 





Сечение А-А' 


Рис. 1. Ориентация осей деформационных смещений и геометрические соотношения, 
определяющие формирование площади 5 срезаемого слоя 
Считается, что жесткость подсистемы заготовки на порядок больше, чем жесткость инструмента, поэтому ее 
деформационными смещениями пренебрегаем. Заметим, что учет упругих деформационных смещений заготовки не 


представляет сложности. Параметры Ш, Й и С можно определить по правилам, изложенным, например, в [11]. 
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Главная проблема при моделировании динамической системы состоит в представлении Р (Г) ={Ё (1), 2, (1) |а 


в координатах состояния. Будем следовать известным взглядам на динамическую связь, формируемую процессом ре- 
зания. Отдельные элементы этих подходов изложены в выше цитированных работах. В целом, они базируются на сле- 
дующих особенностях формирования сил резания: 

Силы возрастают при увеличении площади срезаемого слоя 5; 

— существуют запаздывания сил по отношению к вариациям упругих деформационных смещений [3, 6-9]; 

— они изменяются при варьировании скорости резания [4]; 

— при сближении задней грани инструмента и обработанной части заготовки имеет место непропорциональ- 
но быстрое нарастание сил, действующих на заднюю грань инструмента [6, 10]; 

— существует запаздывание тангенциальных составляющих сил по отношению к их нормальным составляю- 
ЩИМ. 

Кроме этого при моделировании сил используется гипотеза о неизменности ориентации сил в установившем- 
ся состоянии. Она определяется: 

— для главной составляющей силы Е“ угловыми коэффициентами | и Хо; 

— для силы, формируемой в области контакта задней грани с обработанной частью поверхности, — угловы- 


ми коэффициентами у” их’. 


Коэффициенты в основном зависят от переднего угла '/ инструмента (угловые коэффициенты Хи д.) иот 


1 1 ы 
заднего угла (угловые коэффициенты х' и у.) Этим свойствам соответствуют следующие представления сил 


резания в координатах состояния системы: 
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ХО Б® окр 
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[М 
Г О нЕ А 
Ау 
Г, = У =Х,, 
где а, Ь — соответственно толщина и ширина срезаемого слоя, мм; Т,, Т, — постоянные времени, моделирующие 


запаздывание сил по отношению к деформационным смещениям, с, причем, Т,)Т,; р — давление стружки на перед- 


0 ы 
нюю грань инструмента в кг/мм*; Е } — сила резания, действующая на заднюю грань инструмента в состоянии рав- 
новесия, кг; Г, — параметр, определяющий крутизну нарастания силы, действующей на заднюю грань инструмента 


(зависит от заднего угла инструмента 0), мм/с; У — скорость резания; и, ©, — параметры кинетической характе- 


ристики процесса резания, определяющие тангенциальную составляющую силы, действующей на заднюю грань ин- 
струмента, кг-с/мм и с/мм? соответственно. 

Здесь зависимость главной составляющей силы резания от скорости аппроксимируется в соответствии с 
предложениями [4]. Эта составляющая также ориентирована в плоскости с помощью угловых коэффициентов Хх” И 


(2) 
2 


Таким образом, исследованию подлежит следующая система, которая учитывает все отмеченные выше осо- 


бенности динамической модели процесса резания: 


Бир://уезциК.Аопзва.ги 


13 


Механика 


14 


Заковоротный В.Л. и др. Бифуркации стационарных многообразий 





42°, ах ах, 
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2 1 2. ы 
т РК +Й 2 я о Е +< 21 +62 2Х> (3) 


АХ, / 4 АХ АХ 
Е (0 (0) (2) 2 2.3. 
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Главная проблема, рассматриваемая в настоящей статье, связана с анализом притягивающих многообразий, 
которые формируются в окрестности равновесия системы (1), (2). Поэтому вначале рассмотрим свойства равновесия 
— прежде всего, условия потери его устойчивости. 


«р 2 
Свойства равновесия системы. Точка равновесия ( Ху’, Х› ) определяется из условия 4`Х _@Х _(: 
Аг 4 





* * 
+ Р.В) +е.Х. 


(с 
1.1 (4) 


+ Хор-ВХ, че _Х* р.Ь-а ЕО Е о а. (И)?1. 


(С 2 и. 


Система (4) является линейной, и она задает единственную точку равновесия, то есть при изменении парамет- 
ров и внешних воздействий ветвления точек равновесия не наблюдается. Определение ее устойчивости, как известно 


* * 
[26], основано на линеаризованном в окрестности равновесия уравнении в вариациях относительно ( Х г Хх › ) после 
* * <“ 
замены Х\ (К =Х, +х (@ и Х, (1) ЕХ, +х, (1) . При этом учтем, что при анализе устойчивости точки равновесия 


х (1) > 0, х, (1) > 0, ах ()1&=5 0, ах, (1) 14! > 0. Тогда при выполнении этих условий справедливо 


2, 
в +Й; и 0, (5) 
аг 
п, РР ПЕ п. Ей, И си + ре с. 


Т 
где Й; 5 х={х,х,} — вектор ва- 


В. -ХоРЬТ, +ХФЕ® ИТ, в, -ХФа | въ -+хьрЬ с 
— Х2Р Х.2 3 22 Х2 И 1,2 Т Х2Р 20 

риаций упругих деформационных смещений. 

Из (5) получаем характеристический полином линеаризованной в окрестности равновесия системы 


тр” + (1, -ХРЬТ, +? Е® /Т,)р+ (с. +хр-Ь) (Ха )р-+с,, 


ДА(р)= (6) 


(0 р" (0) 2 Хх 
(Й2 — Х2рБТ, + Е ИТ )р+(с,. +Х›р-6) тр” + (й,. — о) р- с, › 
Очевидно, что условия потери устойчивости определяются как свойствами динамической подсистемы ин- 
струмента, так и параметрами и моделями представления сил в координатах состояния системы. В частности, если 


ы я» (2) 
конструктивно обеспечить динамическую систему инструмента ортогональной, то при Й, =) ©. система будет 


устойчивой, если выполнены следующие требования: й,).р.В.Т, — ХЕ ео 5 й. 5) ба 


Ранее показано [11], что недиагональные элементы матрицы С зависят от ориентации эллипса жесткости 
подсистемы инструмента в рассматриваемой плоскости. Этот случай, несмотря на некоторую его условность, подчер- 
кивает, что главные факторы, вызывающие потерю устойчивости, связаны с формированием запаздывающих аргу- 
ментов в преобразовании деформационных смещений в силы резания. Они указывают также на значение угла наклона 
скоростной зависимости сил по мере увеличения скорости резания. 

В общем же случае главный механизм потери устойчивости при резании связан с образованием ускоряющих 
сил, которые формируются симметричной составляющей матрицы скоростных коэффициентов. Суммарная матрица 
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ь 0,С ь 
й; в (5) не симметрична. Следовательно, она представима в виде суммы симметричной и и кососимметричной 


р“ составляющих. Из (5) необходимым условием устойчивости является 
(0.©) _ [1 -Хр-В- 1, + у В” /Ту)] [, › — 0,5(.р-ТЬ + 1” = хо /Т.)] 
: [Я, › —0,5(%›р-Т, + Ха — о /Т. | (й., — ха) 


из которого вытекают все известные условия потери устойчивости. 


У, (7) 


Преобразованию матрицы во из положительно определенной в отрицательно определенную способствует 


увеличение: 
— запаздывающего аргумента Г’; 
— коэффициента <, , определяющего «падающую» характеристику приращения сил при увеличении скоро- 


сти резания. 
В результате система может потерять устойчивость равновесия, так как гироскопические силы, формируемые 


матрицей 
0 Об Ро 
НО 0 


(0,К) _ 
> = 


0,С Г 
не могут стабилизировать равновесие при условии р; (0 [27, 28]. Характеристический полином (6) позволяет в 


параметрическом пространстве проанализировать области устойчивости. Однако нас в данной работе в основном бу- 
дут интересовать условия, при которых система теряет равновесие и в его окрестности формируются различные ста- 
ционарные многообразия. Проанализируем этот вопрос на примере. При этом учтем, что зачастую примеры являются 
более показательными, чем анализ общего уравнения. 
Бифуркации стационарных многообразий в окрестности равновесия. Система (3) моделирует два взаимосвязан- 
ных колебательных контура, которые имеют два источника самовозбуждения. Один связан с запаздыванием сил реза- 
ния по отношению к деформационным смещениям. Второй обусловлен существованием участка вариаций скорости 
(причина его возникновения — так называемый температурно-скоростной фактор), на котором увеличению скорости 
соответствует уменьшение сил резания. Поэтому в системе в зависимости от параметров могут существовать три ста- 
ционарных состояния. Первое — асимптотически устойчивая точка равновесия. Второе — орбитально асимптотиче- 
ски устойчивый предельный цикл. Третье — двумерный инвариантный тор. 

Рассмотрим систему, параметры динамической модели инструмента которой приведены в таблице. Достаточ- 
но малые значения элементов матрицы жесткости выбраны для того, чтобы увеличить чувствительность формируе- 
мых многообразий к изменениям параметров динамической связи в процессе резания. Продольное точение стали 20Х 


осуществляется четырехгранными неперетачиваемыми пластинками из Т16К6б. Основные геометрические параметры 
0 
инструмента: у=6', а=2°, ф=45°. 
Таблица 1 


Параметры динамической модели инструмента 


т, кгх с”/мм В, кгхс/мм с, кг/мм 


0 0,9810 2,0 12 25 108 


Тот 0 30 2,0 т 9 





Частоты системы без взаимного влияния колебательных контуров (Й,, =Й,, =С. =С., =0) и динамиче- 
ской связи, формируемой процессом резания, соответственно равны: 0“) = 1020Гиц, 60), = 332 Гц. Частоты системы 
-, -, 2 2 
с учетом динамической связи, формируемой процессом резания — 60 = 15401 цу, 00 = 840 Гу. Давление стружки на 
переднюю поверхность инструмента — р = (450 —600)кг/ мм’. Коэффициенты ориентации силы резания в плоскости 
авны: Хх, =0,6, х, =0,8, Х® =1, ХО =0,® =0,8 , Хо =0,6. О 
р : Х. =0,6, Хх. =0,8, ХХ =1, ХХ. =0,Х: =0,8 ‚, Х; =0,6. Они удовлетворяют условиям нормировки, 


например (5, )” +(%,)” =1. 
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Чтобы варьировать основные параметры динамической характеристики процесса резания изменялись техно- 
логические режимы: 

— скорость резания в пределах 60-160 м/мин; 

— глубина резания в пределах 0,5-3,0 мм. 

Величина подачи на оборот выбрана неизменной и равна 0,1 мм/об. Все параметры идентифицированы по 


правилам, изложенным в [11]. 
Проанализируем подробно бифуркационную диаграмму в плоскости двух параметров ( Т,- 0), приведенную 


на рис. 2. 
На этой иллюстрации в качестве примера показаны точки 1, 2, 3, 4, для которых на рис. 3 приведены фазовые 


траектории в сечениях фазового пространства. Части а, в, д, ж рис. 3 соответствуют сечению фазового пространства 


Хх 

АЛИ Части 6, г, е, з— сечению Х_—”. Все фазовые траектории сверху вниз соответствуют точкам 1, 2, 3, 4. Также 
ба а 
| 1 


на рис. 2 приведены области О-разбиения для линеаризованного в окрестности равновесия уравнения. Они показаны 


Х 


пунктирными кривыми. На участках А — О — О пунктирные кривые совпадают с границами бифуркационной диа- 
граммы, поэтому они сливаются с границами четырех областей О-разбиения: 

— область, в которой все корни характеристического полинома расположены в отрицательной комплексной 
полуплоскости (А — О — О); 

— области, в которых пара комплексно-сопряженных корней расположена в правой комплексной полуплос- 
кости (А— О—Еи Е — О — О); 

— область, в которой две пары комплексно-сопряженных корней расположены в правой комплексной полу- 


плоскости (Е — О— В. 


Г, мм/с 
50 





45 


40 


35 


30 


0.15 0,2 


о, Кг:с/мм 
Рис. 2. Бифуркационная диаграмма в плоскости (Т., (©): 


А— О — О — область устойчивости равновесия; В — О — С — область формирования двумерного инвариантного тора; 
В — область формирования орбитально асимптотически устойчивого предельного цикла на низкой частоте; 
Р — область формирования орбитально асимптотически устойчивого предельного цикла на высокой частоте 














Как видно, границы областей О-разбиения не совпадают с границами бифуркаций инвариантных многообра- 
зий в параметрическом пространстве. 

Для объяснения этого смещения границы ОЭ-разбиения и линии бифуркации перехода от предельного цикла к 
инвариантному тору дальше можно рассмотреть фазовые траектории, соответствующие точке 2 на рис. 2. 
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Рис. 3. Фазовые траектории системы, соответствующие четырем точкам 
на рис. 2: 1 (а), (6); 2 (в), (г); 3 (9), (е); 4 (ж), 5 (3) 


Временные диаграммы для данного случая приведены на рис. 4. Здесь хорошо видно, что двухчастотный про- 
цесс с течением времени преобразуется в одночастотный, то есть имеет место затягивание высокочастотных колеба- 


ний к низкочастотным. Аналогичная ситуация возникает и в области С — О — Е. Однако здесь низкочастотный про- 


Бир://уезциК.Аопзва.ги 


Механика 


18 


Заковоротный В.Л. и др. Бифуркации стационарных многообразий 





цесс стягивается к высокочастотному. Эти эффекты связаны с нелинейными диссипативными свойствами системы, 


которые учитываются в правой части системы (3). 


0,045 
Х,ь мм 0,04 


0.035 





Рис. 4. Пример преобразования двухчастотного процесса в одночастотный, 
то есть в орбитально асимптотически устойчивый предельный цикл 


Приведем также спектральные характеристики колебаний (рис. 5, а, в, д) для точек 2, 3, 4 на рис. 2. 


Важно подчеркнуть, что при этом меняется расположение ориентации колебаний в плоскости (Х, -Х,)и их 


размах в направлении Х, , в котором инструмент оставляет следы на обрабатываемой поверхности (см. рис. 5, 6, г, е). 


При переходе от низкочастотных автоколебаний к высокочастотным размах Х, соответственно уменьшается с 
0,065 мм до 0,006 мм. При этом изменяется ориентация направления колебаний в плоскости (Х/ -Х,). 


Анализ влияния упруго-диссипативных параметров подсистемы инструмента на формируемые многообразия 
показывает, что качественно бифуркационная диаграмма остается соответствующей рис. 2. Варьирование параметров 
матрицы С приводит к смещениям корней характеристического полинома (6) в комплексной плоскости. 

При этом корни могут сближаться. В таком случае область, в которой формируется двумерный инвариантный 


тор, вырождается. На формирование этой области оказывают существенное влияние недиагональные элементы мат- 


рицы жесткости. При увеличении прежде всего диагональных элементов матрицы Й область устойчивости расширя- 
ется. Таким образом, при изменении параметров упругости и диссипации подсистемы инструмента происходят изме- 


нения многообразий, формируемых в окрестности равновесия системы. 
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Рис. 5. Изменение спектров установившихся колебаний и их ориентация 


в плоскости (Х, —Х,) для точек 2 (а), (6), 3 (в), (г), 4 (9), (е) на рис. 2 


Заключение. Динамическая система резания, формируемая на основе взаимодействия подсистемы инструмента с ди- 
намической связью процесса обработки, является принципиально нелинейной. В окрестности равновесия такой си- 
стемы, как правило, формируются инвариантные многообразия типа предельного цикла или инвариантного тора. Па- 
раметры этих многообразий обусловлены свойствами формируемой динамической связи, зависящей от технологиче- 
ских режимов и от динамических свойств подсистемы инструмента. В свою очередь, известно [6, 29], что, например, 
стойкость инструмента обусловлена параметрами его автоколебаний. При этом в зависимости от частоты имеет место 
оптимальное значение автоколебаний, при которых стойкость максимальна. Кроме этого, в зависимости от ориента- 
ции параметров автоколебаний в пространстве они могут приводить к уменьшению величины шероховатости форми- 
руемой поверхности. Поэтому раскрытие механизмов образования и прогнозирование их параметров открывает путь 
улучшения процесса резания. 
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